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Abstract 



\^ , We construct the Hadamard Green's function by using the eigenfunction, 

which are obtained by solving the wave equation for the massless conformal 

, scalar field on the S" 1-1 of a n-dimensional closed, static universe. We also 

00 . 

, consider the half space case with both the Dirichlet and the Neumann bound- 

ary conditions. Solving of eigenfunction and eigenvalues of the corresponding 



field equation is interesting since the Casimir energy could be calculated an- 
alytically by various methods. 

I. INTRODUCTION 

In 1948, Casimir showed that one consequence of the zero-point quantized field is an 
attractive force between two uncharged, perfectly conducting parallel plates. The principal 
message of this force is that changes in infinite vacuum energy of the quantized field can be 
finite and observable P,|3[|. 

The Casimir energy leads to negative energy between the plates, depends on the distance 
between plates and independent plate well as similar calculation were made for 

sphere [|J, for which the Casimir energy was found to be positive which means that the 
force on the shell is repulsive and depends only on the radius of the shell. The type of 
Casimir energy values depend on manifolds of different topology and geometry ||. The 
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boundary and/or curvature conditions play the important role of the quantized field. 

The Casimir effect has also been calculated in curved background spacetime. In 1975, 
Ford || showed that the vacuum energy of a massless conformal scalar field in an Einstein 
universe by the mode sum method. Ford obtained the renormalized vacuum energy density 
in this case as 

he m 
9 480^' W 

where, Ro is the radius of the universe and the pressure is given by P = |p. 

Thus the energy momentum tensor is in the same form as that for the classical radition. 
In Ford's work on the Casimir effect in an Einstein universe, he showed that the energy is 
associated with the closed spatial topology. Ford used the mode sum method and employed 
an exponential cutoff in his calculations. Later, Dowker and Critchley J7j verified Ford's 
result by using the covariant point-splitting method. Their calculation was not only covariant 
but also cutoff independent. Even though Dowker and Critchley used the global topology 
of the universe to construct the Green's function, their result is also expected to be true 
locally; i.e., the local renormalized energy density and pressure should be given by 

<«> => > = «P5j ■ (" W> = ~ ^ = ~ ^ =) P 4 ' < 2 > 

Three often used regularization techniques of the quantum vacuum energy in a curved 
background spacetimes are the mode sum method 0], the zeta function regularization tech- 



nique |10| and the point-splitting method. One of them, the must powerful method available 
is the point-splitting regularization technique |||9| . This method works with point separated 
functions. In the Hamiltonian, the scalar field <£> (x) behaves like a field operator due to the 
canonical quantization scheme. However, the products of field operators and their covariant 
derivatives are at the same point which causes the vacuum expectation value of the energy 
momentum tensor to the diverge. To eliminate these divergent quantities, one of the field 
operators $ (x) in each product in the Hamiltonian replaced by $(V), where x' is some 
point near x. Thus the separeted function (0 |$ (x) $ (x') \ 0) transforms like the product two 



functions, one at x, the other at x' . The two point function (0 |$ (x) $ (x') \ 0) is called the 
positive Wightmann function and is denoted by 

G {+) (x,x') = (0|$(x)$(x')|0). (3) 

The vacuum expectation value of the anticommutator of the field is called the Hadamard 
Green's function, which is defined by M 

(x, x) = (0 |{$ (x) $ (a/) + $ (a/) $ (x)}| 0) , (4) 
where (x, x') and (x, x') are related by 

G {1) {x,x') = 2ReG {+) {x,x') . (5) 

In this work, we use the above physical motivation for the massless conformal scalar field 
in a n-dimensional closed, static universe. Firstly, we will construct the Hadamard Green's 
function for the case of a n-dimensional static, closed universe, using the separated function 
idea. Secondly, we construct the Green's function directly by using the wave function satis- 
fying the boundary conditions. The half space case (as denoted by Kennedy and Unwin |TTJ) 
is very important in the sense that despite being in curved spacetime, and with a spherical 
boundary. Our results confirm the Green's function for the full n-dimensional closed, static 
universe as an image sum and obtained the necessary Green's function for the n-dimensional 
half space case by locating an image charge in dual region. 

We have organized this paper as follows. In section II. we calculate the mode function 
for the massless conformal scalar field in a n- dimensional closed, static universe. In section 
III. we constructed the Hadamard Green's function by using the eigenfunction, which are 
obtained by solving wave equation for the massless conformal scalar field on 5* n_1 (n > 4, 
n is the dimension of the spacetime). In section IV. we study the Green's function for the 
massless conformal scalar field with a spherical boundary at Xo — § and with the interior 
geometry represented by the closed, static n-dimensional universe (this case is also called 
the half universe). We consider both the Dirichlet and the Neumann boundary conditions. 
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II. THE MODE FUNCTION 



We shall use hyperspherical polar coordinates (x,9,4>) = (x, #i,#2,#3, 1^0) ! I 1 — 

1,2, 3, , n — 3. Where n is the spacetime dimension. The line element in the time orthog- 
onal form 

ds 2 = dt 2 - R 2 dl 2 , (6) 

where Ro is the radius of the universe and is a constant. In a static-closed model the 
part dl 2 follows by imagining that the space sections of fixed world time t are embedded in 
a n-dimensional space with coordinates x±, x 2 , x 3 , x n _i, w and line element dl 2 = dx\ + 

dx\ + + dx 2 n _ x + dw 2 . Our n — 1 dimensional space is the surface of the sphere x\ + x\ + 

+ x 2 n _ x + w 2 = Rq at fixed R . Spherical coordinates for this space are 

X\ = R sin x sin 9i sin 9 2 sin 9 n _ 3 cos 0, (7) 

x 2 = Ro sin x sin 9 X sin 9 2 sin 9 n ^ 3 sin 0, 

x 3 = R sin x sin 9 X sin 9 2 cos n _3, 

x n _i = R sin x cos 9i , 
w = Ro cosx- 

Where x e [0,tt], 9^e[0,n] /x = 1, 2, 3, n - 3, and <p E [0, 2tt]. 

Thus, the general metric for a static, closed n-dimensional spacetimes has the form 

ds 2 = dt 2 - R\\dx 2 + sin 2 xd9 2 + sin 2 x sin 2 9 x d9\ + 

+ sin 2 x sin 2 9\ sin 2 # 2 sin 2 9 n -±d9 2 n _ z 

+ sin 2 x sin 2 9\ sin 2 6*2 sin 2 # n _4 sin 2 ^ n _ 3 rf0 2 ] . (8) 

We shall consider conformal massless scalar field $ (x) on this curved background geom- 
etry. It satisfies the field equation 
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□$ (x) + £ (n) i?$ (x) = 0. 



(9) 



(We follow the sign convention in ||) Where £ (n) = | fc-i) is t ne dimensional parameter 
determining the coupling between $ (x) and the Ricci scalar R. The Ricci scalar for the 
metric Eq.(8) is given by 



R = — {n - 2) ( n - 1) for n > 2. 

Rn 



(10) 



We solve Eq. (9) using the method of separation of variables. In general, for n > 4 we 
can write 



$ (x) = c 0n e~ lM X ( x ) X 1 (6 1 ) X 2 (0 2 ) X n _ 4 (0„- 4 ) (0„- 

with o> = — . 



(11) 



Where Y\ m is the spherical harmonic and co n is appropriate normalization constant, and 



(N 2 =)/ = 0,1,2,3, ,and m = -I, -I + 1, 0, 1, (12) 



We now calculate the mode function for n > 4. The X (%) and X^ (6^) (/i 
1,2,3, ,n — 4 ) satisfy the following differential equations 

sm xdx V «X 



+ 



iV 2 - 



("~2) 
4 



sin 2 x - iV„- 2 (iVn-2 + n - 3) 



X( X )=0 ,and (13) 



1 L( sin n-^2 9 d 

sin n_At 6,i d9u \ " W '' 



+ [iV n _ M _ 1 (X n _ M _! + n - /i - 2) sin 2 # M - N n ^ 2 (N„_ M _ 2 + n - // - 3)] X„ (6> M ) = 



(14) 
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with fj, = 1,2,3,4, ,n- 4. 

Where ( N 2 =)l, N 3j N 4: N 5 , N n _ 3 , N n _ 2 and N = R u are the separation parame- 
ters. 

Notice that we still have not calculated the regularity of the solution in our intervals 

at # M e[0,7r] where [i = 1,2,3, ,n — 4; hence Nj (j — 3, 4, 5, n — 2) and TV = Rquj 

remain continuous parameters while I and m take the values (N 2 —)l = 0, 1, 2, 3, 4, 5, 6, , 

m — — I, — I + 1, 0, 1, ,1 — 1,1 . 

Making the following substitutions for each one in Eqs. (13) and (14) 

X ( X ) oc sin^- 2 X ( X ) Ci (cos X ) , and 



X^O^) oc sin^-"- 2 ^ C M (cos0 M ) 
with fi — 1,2, 3, 4, , ?7,-4. 



(15) 



We obtain 



(1 - x 2 ) -^C (x) - x (2iV n „ 2 + n - 1) -^-C (x) 

ax z ax 



+ 



- N n „ 2 (N n „ 2 +n-2) 



C(x) = , 



(16) 



where x = cos X . 

d 2 d 
i 1 ~ x l) 7^ C » M ~ x » ( 2iV «^-2 + n-n-l) ^t^m M 



dx 2 

H 1 



+ [W„-„-i (A^n-M-i + n - /i - 2) - iV n _ M _ 2 (iV„- M _ 2 + n - /i - 2)] C M (x M ) = 0, (17) 



where x^ = cos 0^ , fj, — 1,2, 3, 4, , n — 4 . 



The differential equation Eq. (17) can be solved by the method of Frobenius. All solution 
sets are also divergent at the end points of our interval. To obtain regular solutions in the 
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entire interval e [—1, 1] (/i = 1, 2, 3, 4, , n — 4) we terminate the infinite series after 

a finite number of terms by restricting N 3 , N±, , N n _ 2 to have integer values. 

7V n _ 2 = 0,l,2,3,4, , 

Nn-n-i = 0, 1, 2, 3, 4, , Nn-n with ft = 2, 3, 4, n - 4, 

(jV 2 =)Z = 0,l,2,3,4, ,7V 3 

m = _/ > _/ + i j ; 0,1, ,/ . (18) 

The polynomial solution obtained this way for the equation (17) is the Gegenbauer 
polynomials 



X M (^)oc sin"— ^C^:;V_ m _ 2 (cos^) with ^ = 1,2, 3,4, ..,n- 4. (19) 

We return to Eq. (16). This differential equation can be solved by the method of 

Frobenius. Hence, we express the function C (x) by a power series as 

00 

C(x) =J2 a k x k+a . (20) 

fc=0 

This leads us to the following recursion relation: 

^2 



«fc+2 



k + iV n _ 2 + 



— N 2 

. (21) 



a k (k + l)(k + 2) 

The above recursion relation gives us the following series solution for C (x) : 

00 00 

C{x) = Y, a 2k x 2k + «2*+i x 2k+1 . (22) 

k=0 k=0 

Using the Recurrence relation this is in the form 

C(x) = a C 1 (x)+a 1 C 2 (x) , (23) 

Where C\ (x) and C 2 (x) are linearly independent. Convergence of the series at the end 

00 

points of our interval can be checked eaisly by using the Raabe test which says that if 

k=0 

uu is a series of positive terms and if the lim n ( — l) = A the series is divergent for 



A < 1 , convergent for A > 1 and the test fails for A = 1 [T5]. At the end points we apply 
this test to the series with even powers i.e. C\ (x) . We write 

a 2k (2k + 1) (2k + 2) 



2A: + 7V„_ 2 



("~2) 
2 



(24) 



iV 2 



one obtains A 



N, 



n-2 



(n-2) 



(with N n _ 2 > and n > 4) A < 1 . Hence the 



infinite series C\ (x) diverges at the end points. One may similarly check that C 2 (x) is also 
divergent at the end points of the interval. To obtain a regular solution in the entire interval 
x £ [—1,1] we terminate the infinite series after a finite number of terms by restricting N 
to have integer values given as 

(n-2) 



N(k,N n „ 2 ) = k + N n _ 2 + 



(25) 



where k = 0, 1, 2, 3, 4, and N n _ 2 = 0, 1, 2, 3, with u 



Ho 



. Redefining a 



new index /c' as k + iV n _2 we write N = k' + 2 with u>k = jf and fc' = 0, 1, 2, 3, 



and iV n _2 = 0, 1,2, 3, Dropping primes one obtains the eigenfrequencies as 



k + 



[n-2) 



with k = 0,1,2,3,4, 



(26) 



Values that N (k, N n _ 2 ) could take are presented in table I.. 
For fixed k, one obtains the degeneracy of each eigenfrequencies 

k N n - 2 N 6 N 5 N4. N 3 

* E E EE EE «• '• 

7V„_ 2 =0 7V n _ 3 =0 AT 5 =o 7V 4 =0 V 3 =0 Z=0 



(27) 



The polynomial solutions obtained this way are the well known Gegenbauer polyno- 
mials, hence we can write G\ (x) in terms of the Gegenbauer polynomials as C\ (x) = 

N 1 ("-2) 

Cfc-^L 2 ( cos *)- 

Thus, the mode solution for n > 4 becomes 



$ a (s) = c „ e"*** sin^-x CfcZI, 2 ( cos x) 



n— 4 



n[ 



AT 1 n — fi — 2 
JVn-n-2H 1 



iV, 



n — ^ — 1 -^n — fi — 2 



(cos#„ 



m ^"n— 3, 



(28) 



with 





1 

Rq 


[ k+ 2 


2)" 


, where k = 


0,1,2,3,4, 


N n . 2 = 


0,1,2,3,4,.... 


..,k 




N n -(fi+i) = 


0,1,2,3,4,.... 




..,N n n with 


A* = 2, 3, 4 


(N 2 =)l = 


0,1,2,3,4,.... 




....,N 3 and 




m = 


-1,- 


-/ + !,.... 


•,0,1,2, ,/ , 





(29) 



and con is the normalization constant. 

To evaluate the normalization constant we write the scalar product 

($ Al $ A2 ) =-iJ [$ Al (x) d t $* 2 (x) - d t $ Al (x) $* 2 (*)] d"" 1 * 
= 5\ 1 \ 2 . 



(30) 



$ A (A stands for k, iV n _ 2 , AT n _ 3 , , A^, ./V3, (AT 2 =)/ and m values which are given in Eq. 

given Eq. (28) and, from a complete orthonormal set with respect to the product Eq. (30). 
Using Eq. (30) we obtain co„ as 



("-4) 



COn — 







n-2) 



(k-N n _ 2 )\ 



n— 4 



n * 



. (n— tt— 3) 

-„-2+— 5— r Nn _^_ 2 + 



(iV„- M - 



yrr (k + N n _ 2 +n-2) 
(n-n-2)\ 



vrr (A^„_ M _i + A^„_ M _2 + n - n - 2) 



(31) 



The general solution of the field equation can be written as sum over these modes in the 
form 



$ (X) = J2 $ A (X) + 4$ A (x)) 



(32) 



Where, a A is an annihilation operator and a A 's a creation operator for quanta in mode 
A, while $ A (x) is the solution of the wave equation given by Eq. (28). 
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III. THE GREEN'S FUNCTION 



The Green's function we aim to calculate is defined by 



(0 | $ (x) $ (x')\ 0) = $ a (x) $1 (x') 



(33) 



where $a (x) are the mode function. 

Using the mode solutions given in equation one obtains 



(x,x') 



irK 



—y 

m-2 
fc=0 



sin Afn-2 xsin JV„_2 x / Q 



N n - 2 =0 ^ ^ 



(n-2) 



r (N n _2 + k + n-2) 

n— 4 I Nn-^-i . . 

At=l [Ar n _ M _ 2 =0 ^ ^ 

vrr {N n ^ + A^_ 2 + n - n - 2) 

(n-M-2) 



cosx) <?jfe-jvL a 2 (cosx') 



(n- fi-2) 



k-N n -2 
2 



(n— /t— 2) 



sin^-- 2 ^ sin^--^; Cg:^^ fcosfl,,) G " ^ ^ 



A r „- M _ 2 + i 



(cos * y 



3) 



(34) 



m=—l 



Using the addition formula for the spherical harmonics |L2]], which is given as 

/ 

— Pi fcOS "Yi) = 

47T 



^±l)p,(cos 7l ) = J2Y lm (£ n _ 3 ,0) ^(C 3 >0 , 



m=— I 



where cos ji = cos 6 n ^ 3 cos 6^_ 3 + sin 6> n _ 3 sin #^_ 3 cos (</> — 0') 



(35) 



And, using the n — 3 times addition theorem for the Gegenbauer polynomials |i~3" | 



r<P( i I • n ■ / \ r (2p - 1) ^ 2 2 (z/-r^)! (2?7 + 2p- 1) 

Gf; cos w cos (/? + sm 17 sin w cos w = tt- > 2 ' T p + n : : 

{r(p)} 2 £^ r(z7 + ?7 + 2p) 



??=0 



sin^inV C^(costf) C^(cosp) 5 (cos^). 



(36) 
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We obtain 



G (+> (*• = 2535^ E «"«™ ^ (cos 7 „_ 2 ) , (37) 

where 



a =2 o 2l22 ^ 4 {r(i)} 2 {r(l)nr(2)} 2 {r(^)y 

r (i) r (2) r (3) r (4) r (n - 3) ' 



(38) 



and 



cos 7„_ 2 = cos x cos x' + sin \ sin x' cos 7„_ 3 , 
cos 7„_^_ 2 = cos 6> M cos 6^ + sin 6^ sin 6^ cos 7„-^-3 
with //= 1,2,3,4, ,n — 4 , and 

cos 71 = cos 0„_ 3 cos 0' _ 3 + sin 6> n _ 3 sin _ 3 cos (</> — 0') . (39) 

Substituting At — ie for with the understanding that we will let e — * at the end we get 
< 1 . Now, we could use the generating function for the Gegenbauer polynomials, 
one obtains 



e R o 



£(+) U x >) = ^ . (40) 

/ At \ 2 

[cos - COS 7„_ 2 J 



Where 7„_2 is the angle between the vectors r and r' in the direction 

(X, 0i, fc, 03, ,0n-3,0) and (x',0i,0 2 ,03, A-3,0') ■ The geodesic distance on S"' 1 

is denoted by As n _ 2 (g, <?') = Roln-2 where q = q(x,0i,9 2 ,0 3 , , n _ 3 , 0) ; g' and g G S"^ 1 

and x = x (t,q) . The spacelike separation As n _ 2 (g, q') is given in terms of the coordinates 
(X, 0i, 02, 03, , 0„-s, 0) and (x', 0i, 2 , 3 , , 0U, 0') as 

(cos 7„_ 2 =) cos— — = cos x cos x + sin x sin x cos 7„_ 3 , (41) 

Kq 

where cos 7„_ 3 is given in Eq. (39). 

Thus, we have the following the two point function 
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(° l $ (*) $ Ml 0) = _ 2 - ^ . (42) 

" ^cos#-cos A - 



Notice that the Hadamard Green's function (x, a/) is related to (x, x') through 
the relation 

G m (x,x') =2Re (0|$(x) $(x')|0) , 

= 2Re G {+) (x, x') . (43) 



IV. THE GREEN'S FUNCTION FOR HALF UNIVERSE 



We will now investigate the massless conformal scalar field with a spherical boundary 
located at xo — f an d with interior geometry represented by the static closed metric. The 
metric inside the boundary may be expressed as 

ds 2 = dt 2 - Rl[dx 2 + sin 2 \dQ\ + sin 2 x sin 2 9 x d9 2 2 + 

+ sin 2 x sin 2 9\ sin 2 9 2 sin 2 9 n -$_d9 2 n _ z 

+ sin 2 x sin 2 #i sin 2 9 2 sin 2 # n _ 4 sin 2 9 n ^d(j) 2 ] , (44) 



where x e [0, f] , # M e[0,vr], /i = 1, 2, 3, 4, 5, n - 3 and e [0, 2tt] . 
We consider a massless conformal scalar field on this curved background geometry where 
the wave equation is given by 

n$(x)+£(n)R$(x) = 0. (45) 

Where £ (n) = | ^"~ 2 | is the conformal coupling constant and i? is the Ricci scalar: 
R = \ {n - 2) (n - 1) for n > 4 (n is the dimension of spacetime). 
Solution of Eq. (45) could be found easily 
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viV„-2 + 



*> x (x) =c£e-** sm N — x C k _ Nn _ 2 

n—A 

n 



("-2) 



(cosx) 



AT„_„._ 2 a ^n-^-l-r 2 



(re-jj-2) 

iV„_ u _ 2 (cos^) 



(#71-3,0) , 



(46) 



where is the normalization constant, and - ,, 

We note that the boundary condition at X = Xo nas n °t imposed on Eq. (46) yet. 
Hence, N still remains a contiuous parameter, while N n _ 2 , N n _^ +1 ^ (fj, — 2, 3, 4, n — 4) , 
I and m take the values. 

7V n _ 2 = 0,l,2,3,4,5, 

iV„- (M+ i) = 0,1,2,3,4, ,iV n _ M , with fj, = 2, 3, 4, n - 4 , 

(JV 2 =)/ = 0,1,2,3,4, iV 3 and 

77i= — Z, — Z h- 1, ,0,1,2, I 



We now impose the Dirichlet boundary condition for our problem i.e. 

7T 



$ A (x) 



X=X0 



at xo 



(47) 



(48) 



Using the wave function given Eq. (46) we immediately get integer values for N, which 
could be tabulated as in table II.. The new eigenvalues ujk = and the degenaracy gu 
could be written as 

1 / (n-2)\ 
UJk= R~ { k + ~2j ' fc = ' 1 ' 2 ' 3 ' 4 '-- 

k — 1, N n _ 2 = 
k = 2, iV n _ 2 = 1 
k — 3, N n _ 2 = 0,2 
fc = 4, AT n _ 2 = 1,3 
fc = 5, iV n _ 2 = 0,2,4 



7V 5 N 4 N 3 

E E E E E • 1 

JV„-3=0 iV„_ 4 =0 JV 4 =0 JV 3 =0J=0 



(49) 



The new eigenfunctions Eq. (46) are still the Gegenbauer polnomials, but the only ones 
that satisfy the boundary condition $ A (x) \ x=xo = are picked: hence table II. could be 
formed by just taking the second, fourth, etc. diagonal elements of the table I.. 
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Now we impose the Neumann boundary condition, which is given by 

7T 



n a V a $ A (x) = at xo 



2 ' 



(50) 



where n a is the outward normal to the boundary. Hence for the Neumann boundary 
condition 



(x) \ X=X0 = at xo = | • 



(51) 



We obtain the N values presented in table III. The new eigenvalues Uk and the degeneracy 
could be written as 

Uk = li ~( k + ^~1T^) ' A; = ' 1 ' 2 ' 3 ' 4 -- 
k — 0, N n _ 2 = 
k = 1, iV n _ 2 = 1 
fc = 2, AT n _ 2 = 0,2 
k = 3, N n _ 2 = 1,3 
fc = 4, iV n _ 2 = 0,2,4 



9k = < 



iv, 



Ti-2 Nn-3 JV 5 7V 4 JV 3 

E E EEEP'+D 

JV„-3=0 iV„_ 4 =0 AT 4 =0Af 3 =0 1=0 



(52) 



Notice that the N values in table III. could be obtained from the full universe case in 
table I. by taking the first, third, etc. diagonal elements. Once again, the eigenf unctions are 
the Gegenbauer polnomials, but only the ones that satisfy the Neumann boundary condition 
are allowed. 

The Green's function that we aim to calculate is defined by 



(0\$(x)$(x')\0) = J2 $ a(^') , 



(53) 



where $a (x) is the solution of the wave equation Eq. (45) with the appropriate boundary 
conditions Eqs. (48) and (51). Using the wave function given in Eq. (46), one obtains 
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(k-N n _ 2 )\ 



AT . (n— 2) * r . (n— 2) 

sin^xsin^V C^:! 9 -^- (cosx) 



T (iV n _ 2 + fc + n - 2) 

n— 4 [ -^n-n-i /- / 

m=i [at„_ m _ 2 =o ^ ^ 

Trr (W n _„_i + A^„_ M _ 2 + n - /x - 2) 



cosx 



(n - /i - 2) 



{n-fi — 2) 



at . (n — /Lt— 2) AT 

A* AT„_ M _l-jV n _ M _ 2 l C0St VJ jV n _ M _l-AT„_ M _ 2 l C0S V 



m=—l 



X(N n ^ 2 ,k) stands the new eigenvalues for the Dirichlet and Neumann boundary condi- 
tions are presented in tables II. and III.. 

We define the Green's function for the Dirichlet and Neumann boundary conditions as 



(54) 



(0\$ D (x)$ D (a/)\0) = D D (x,x') , and 
{0\$ N (x)$ N (x>)\0) =D N (x,x r ) . 



(55) 



We easily see that 



D D (x, x') + D N (x, x') = 2 D (x, x') . 



(56) 



Where D (x, x') = (0 |$ (x) $ (x') \ 0) is the Green's function for the n-dimensional uni- 
verse which given in Eq. (40). Since the sum of the Dirichlet and Neumann eigenvalues 
for the half universe which is depicted in tables II. and III. is equal to full n-dimensional 
universe eigenvalues given in table I.. The factor of 2 comes from the fact that the mode 
functions used in D (x, x') are normalized with respect to the n-dimensional universe. 

To evaluate Do (x, x') and (x, x') explicitly we find it convenient to write them as 



D D (x, x') = D (x, x') — — [D N (x, x') — D D (x, x')\ , and 
D N (x, x') = D (x, x ) + - [D N (x, x') - D D (x, x')\ . 



(57) 
(58) 
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The second term on the right hand side of Eqs. (57) and (58) could be written as 
D B (x, x') = \ [D N (x, x') — D D (x, x')\ . A closed expression for the D (x, x') is already- 
known which given in Eq. (40) and could be written as 

(G (+ Hx,x') =) D(x,x>) = a J - l - . (59) 

47F ^0 ( At As„_ 2 \ 2 

^ cos __ cos __j 

Where At = t - t' , and 

As 2 

cos — — = cos x cos x' + sm X sin x' cos 7«-3 • (60) 
ito 

Here cos7 n _ 3 and a n are given by Eq. (38) and Eq. (39), respectively. 
Using the information given in tables II. and III., definitions of Dn {x, x') and Dp (x, x') 
, and addition theorem for the Gegenbauer polynomials one could write as 

a n T (n — 3) 



D B (x,x') = 



(n-2) (n-2) 

C 2 ( cos x) C 2 (cosx') C 2 (cos7 n _ 3 
+e }T^) sin x sin x' 

n n ("~3) 

C 2 (cosx) C 2 (cosx') c \ 2 (cos7„_ 3 ) 

|p ( r J-l\\ 2 (n-2) (n-2) (n-3) 

- r(n-i) (n ~ 3) Cl 2 (cosx) Cl 2 (cosx ' } c ° 2 (cos7 "- 3) 

+e- 1 ^" | r(n) (n ~ 3) C ^~^^ 

(n-2) (n-3) 

C 2 2 (cosx') C 2 (cos7„_ 3 ) 

2 2 {r(f)} 2 i! 



r (n + 1 ) 



(n — 1) sin x sin x' 



n n (n — 3) 

(cosx) Q 2 (cosx') C 1 2 (cos7„_ 3 ) 

2 4 {r(^)) 2 0' (n + 2) 

+ — w o \ ( n + 1 ) sin2 X sin 2 x'C 2 (cos x) 
I (n + 2 j 

(»+2) ("~3) ^ 

C 2 (cosx') C, 2 (cos 7 n- 3 ) + ] • (61) 



Using the symmetry properties of Gegenbauer polynomials 
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C k-N n . a (-cosx) = (-1) C 



h-N n - 2 



cos/) , 



(62) 



we could now rewrite Eq. (61) as 



D B (x,x') = 



a n T (n — 3) 



oo fc 



r(7v n _ 2 + ^ 



(2iY" n _ 2 + n-3) sirA- 2 x sin^™- 2 x' 



T (iV n _ 2 + fc + (ra - 2)) 



C fc -^_ 2 2 ( cos ^) C k-N n . 2 2 (-cosxj C^ n 2 _ 2 (cos 7n _ 3 ) • 



(63) 



Using the addition theorem for the Gegenbauer polynomials, one obtains 



(n-2) 



47r n o fc=0 



where 



cos 7„_ 2 = - cos x cos x + sin x sm x cos 7„„ 3 . 



e «o 



(64) 



(65) 



< 1 . Now the sum in Eq. (64) could 



Letting At — > At — ie where e — > , then 
be evaluated using the generating function of Gegenbauer polynomials, one obtains 

OL n 1 



D B (x, x') 



cos __ cos _ 



("-2) ' 
2 



where 



(66) 



(cos 7„_ 2 =) cos A p n 2 = - cos x cos x + shi x sin x' cos 7„_ 3 , 
cos 7„_ M _ 2 = cos 9^ cos ^ + sin M sin ^ cos 7 n _ M _ 3 
with /x = 2, 3, 4, n — 4 , and 

cos 71 = cos n _ 3 cos 9' n _ 3 + sin #„_ 3 sin 6' n _ 3 cos (^ — </>') . 
We could now write the complete Green's functions Dm (x,x') and Dr> (x,x') as 



D N (x, x') 



D D (x, x') 



a r 



cos ^ — COS 



(n-2) 
As»_ 2 \ 2 
i?0 



+ 



("~2) 

At As„_n \ 2 

cos __ cos __ 



(67) 



47r n - 2 i?™- 2 



cos ^ — COS 



As„_ 2 \ 
Ro J 



("-2) 
2 



cos f -cos ^ 



("-2) 
2 



, (68) 



(69) 



17 



Where cos A ^ o ~ 2 and cos A ^ o ~ 2 are given by Eq. (39) and Eq. (67), respectively. 

We also note that these Green's functions are naturally identical to the Green's functions 
that one could get by using the image method, which always works on the double manifold 
defined by HHED1 

M U dM U M* . (70) 

M is the physical space inside the boundary, dM is the boundary and M* is the dual 
space obtained by reflecting the physical space about the boundary. The field is confined to 
the region defined by M and satisfies the boundary conditions Eqs. (48) and (50). To satisfy 
the boundary condition on dM it is sufficient to locate an image charge in the unphysical 
dual region M*. The Green's function on M U dM is then given by 

G (x, x') = D (x, x') ± D (x, X*) , (71) 

where — (+) refers to the Dirichlet (Neumann) boundary condition, D (x, x') is the 
Green's function for the double manifold and x' is the image of x'. The key element in 
the success of image method is that the solutions of Eq. (45) found in the double manifold 
will have distinct parity (even or odd) under transformation x — > x. Notice that in Eq. (68) 
and Eq. (69) the first term is the Green's function for the double manifold, which in this 
case is the full S 1 " -1 , and the (cos7^T 2 =) cos A ^'~ 2 in the second term could be rewritten 
as 

(cos t~T 2 =) cos A p n ~ 2 = - cos x cos x' + sin x sin x' cos 7„_ 3 , (72) 
which could be obtained from cos A ^~ 2 by substituting the image of x' i.e. by replacing 



X' ~> 7T - X'- 



V. SUMMARY 

We constructed the Green's function by using the eigenfunctions, which are obtained by 
solving the wave equation for a conformal scalar field in a n-dimensional closed, static uni- 
verse and with the appropriate boundary conditions (Dirichlet and Neumann). Even though 
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we use the global topology of the n- dimensional universe to construct the Green's function, 
our result is locally true. Our result is interesting, to understand the correspondence be- 
tween the sign of the Casimir energy and different manifold structure in closed topology of 
the universe. 

Our result is an agreement with the results obtained by the method of images. The- 
Green's function is calculated for an isolated sphere, where the interior geometry is given 
by Eq. (44) (half space universe). The global topology of the universe as well as the outside 
geometry are not important and have no effect on our result. This very important for appli- 
cations of the Casimir effect to the bag models. Currently we are investigating the Casimir 
effect for a massless conformal scalar field in a n-dimensional closed, static universe and a 
half space with Dirichlet and Neumann boundary conditions. 
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Table I. : N (k, N n _ 2 ) values for the full space case 



N n - 2 \ k 
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n 
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n+4 
2 
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Table II.: N (k, N n _ 2 ) values for the Dirichlet boundary case 



N n .2 \ k 
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Table III.: N (k, iV n _ 2 ) values for the Neumann boundary case. 



N n -2 \ k 
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n-2 
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ra+2 
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